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The Role of Rotor Impedance in the 
Vlbratioa Analysis of Rotorcraft 

Part IV of Final Report under 
Contract NAS2-7613 

As an alternative to a total system vibration analysis with its 
lack of visibility of critical dynamic parameters » one often first 
determines the dynamic rotor forces and moments acting on a rigid 
support, and one then excites the flexible airframe-less-rotor with 
these inputs* In an improved method which still retains the advantage 
of separate treatment of rotor and airframe, the rotor impedance is 
used to correct the input to the airframe* This improved method is 
illustrated for a strongly idealized case of vertical excitation and 
then for rolling and pitching moment excitation of a four bladed 
hingeless rotor on an up-focussing flexible mount* Contrary to the 
usual approach that represents aeroelastic blade motions by a series 
of normal blade modes in vacuum, the aeroelastic rotor impedances are 
computed directly with a finite blade element method that includes 
aerodynamics. The rotor impedance matrix for 3 or more blades is 
determined from the root moment impedance for a single blade by a 
simple multiblade transformation rule* Force and moment amplitudes 
transferred from the rotor to the support are found to be critically 
dependent on the support dynamics. Thus the rotor impedance is shown 
to be an essential item in the vibration analysis* 
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Preface t Part IV of Final Report under Contract NAS2*7613 

Between July 1, 1977 and June 30, 1978 work continued on various 
items covered by the research goals of Contract NAS2-7613. This 
contract that became effective July 1, 1973 was extended several times 
without increase in funding. During FY 1978 it was first extended to 
31 December, 1977, and then to 30 June, 1978. No salary contributions 
were made during FY 1978* The small remaining amount of unexpended 
funds was used only to cover some expenses for typing, computing, 
reproduction, and travel. 

The work performed during FY 1978 covered 3 ^reas: 

1. An extension of Part III of the Final Report under Contract 
NAS2-7613 to the problem of vibration prediction, which is a 
contribution to research goal (a) of the work statement. The 
result of this work is presented as Part IV of the Final Report 
under Contract NAS2-7613. In order to make Part IV self con- 
tained, the finite element method developed in Part III of the 
Final Report is outlined in condensed fom. A paper essentially 
identical with Part IV has been accepted for presentation at the 
4th European Rotorcraft and Powered Lift Aircraft Forum at 
Stress, September 1978. 

2. A second extension of Part III of the Final Report under Contract 
NAS2-7613 to the problem of reducing the order of the charac- 
teristic equation for the coupled rotor/airframe system. This 
proble:n had been only briefly discussed in Part III though it is 
of great practical significance. Work on the second extension of 
Part III is on-going and will probably result in a future publi- 
cation. 

3. An extension of Part I of the Final Report under Contract 
NAS2-7613 to the problem of unsteady wake parameter identifica- 
tion from dynamic model tests, using the hub stirring technique 
rather than the cyclic pitch stirring technique previously 
developed. The hub stirring rotor model has been built and 
operated but needs instrumentation and calibration. If con- 
tinued, this work will also result in a future publication. 

No farther reports unoer Contract NAS2-7613 will be forthcoming. 



iii 


Table of Contents 

Page 

Nomenclature 1 

Introduction 3 

Rotor /Airframe Impedance Matching 4 

Rotor and Support Representation by 3 Masses and 2 Springs 6 

Vertical Force Excitation 7 

Rolling and Pitching Impedance Analysis for a Hingeless Rotor 9 

Blade Element Transfer Matrix 11 

Blade Root Moment Impedance 13 

Transformation to Multiblade Impedance 15 

Coupled Rotor/Support Analysis 17 

Computational Limitations 18 

Coupled Rotor and Support Modes 20 

Hingeless Rotor Moment Excitation 22 

Conclusion 23 

References 24 

List of Figures and Figures 26 to 33 



Nomenclature 


El 

[El 

[F] 

V 
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M 
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N 

P.Q 

R 
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T 

X 

a 

b 

c 

c. 

X 

f/f. 
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EI/EIq, flap-bending stiffness over root value 

transfer matrix across massless blade element 

transfer matrix across point mass 

blade moment over support moment of inertia 

3 2 

blade flap-bending moment, unit R Q 

rotor moment on flexible over moment on rigid support 

number of finite elements per blade 

polynomials of X 

rotor and support impedances 

rotor radius 

2 2 

blade shear force, unit p^R U 

^ 2 2 

blade tension force, unit PqR ^2 
state vector 

blade airfoil lift slope 
number of blades in rotor 
blade chord, unit R 

aerodynamic damping coefficient at 1th station 
rotor force on flexible over force on rigid support 
aerodynamic coefficient for blade pitch angle 
length of blade element, unit R 
one half rotor vertical mass 
support mass 

point mass at station 1, unit p R 

— 1/2 ° 

(T/q El ) blade tension force parameter 

4 2 

EIq/Pq R n blade root flap-bending stiffness 
speed parameter 

distance of blade station i from rotor center, unit R 
time, unit 1/Q 

blade deflection, positive up, unit R 

a single blade variable 

blade pitch angle, positive up 

[E] [F] , transfer matrix 

azimuth angle of kth blade 
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rotor angular speed 
hub tilt angle 

3pac, Lock number for uniform blade 

real part of complex valued frequency, unit Q 

; + lb), complex valued frequency, unit $2 

~ 2 
air density, unit Pq/R 

blade mass per unit length at blade root 

excitation circular frequency, unit Q 

coupled natural frequency, unit (2 

rotor uncoupled natural frequency, unit Q 

support uncoupled natural frequency, unit fi (also (o^, 

support uncoupled natural frequency producing = w 




u)g/o), frequency ratios 


Subscripts 

I, II aultiblade coordinates, forward and left respectively 
i outer end of blade element, beginning with 1 at blade tip 

12, 13 elements of a matrix or determinant 
I, R imaginary, real parts of a polynomial 

M, a, 0 polynomials multiplied by ® 


Superscripts 


time differentiation 
length differentiation 
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Introductlon 


Recently, a great deal of attention has been paid to the problem 
of rotorcraft vibrations. It Is generally accepted that rotorcraft 
vibrations should be alleviated because they reduce crew p'.oflclency, 
cause passenger discomfort, produce equipment deterioration and struc- 
tural fatigue, and Increase maintenance. Fig. 1 taken from reference 1 
shows the cockpit vibration exposure limits at cruise speed In terms 
of hours of maintained crew proficiency vs. frequency In Hz according 
to guidelines In reference 2. The Indicated cruise vibration level 
for UTTAS and AAH shows substantial improvement over their predecessors , 
Che UH-1 and the AH-1. Part of the Improvement Is caused by the 
Increase in vibration frequency when going from a 2 bladed to a 4 
bladed rotor, since humans are less affected by vibratory accelera- 
tions at frequencies above 10 Hz. 

Presently vibration prediction methods are not reliable. There 
is hardly a rotorcraft prototype which did not exhibit during Initial 
flight testing excessive vibrations. Their reduction to specification 
level usually requires an intensive, costly, time and payload consuming 
effort. The vibration level is often sensitive to small variations in 
dynamic parameters so that two rotorcraft from the same production run 
can have substantial differences in vibration level. There are a 
variety of vibration sources. Most important is usually main rotor 
excitation with Che blade passage frequency which will be the topic of 
this paper. Other sources of vibrations are higher rotor harmonics 
particularly with twice blade passage frequency, rotor wake excitation 
of Che empenage, mass, aero, and damper unbalances, rotor self-excitation 
leading to limit cycles, and parametric rotor excitation leading to 
fractional harmonica. Each of these vibration sources requires dif- 
ferent measures of alleviation. For main rotor excitation of the 
airframe the most important means of vibration control is the proper 
selection of the frequency spectrum for rotor and airframe. Other 
means are passive blade pitch control, active higher harmonic blade 
pitch control (not as yet flight tested), rotating system absorbers 
(pendulum, bifllar), fixed system absorbers, preferably self-tuning. 
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and the various means of passive and active rotor isolation. The 
vibration prediction method to be discussed here is applicable to all 
of these dynamic configurations, since they all require a fully coupled 
rotor/airframe analysis. 


Rotor/Airframe Impedance Matching 

Coupled rotor/alrframe vibrations can be and have been treated by 
total system analysis, see for example references 3 to 5. The large 
computer programs needed for such an analysis can easily obscure the 
Influence of crucial dynamic parameters on the vibration level. More 
visibility for such parameters is achieved by performing separate 
vibration analyses (or testing) for the rotor and for the airframe and 
by then matching dynamic forces, moments and deflections at the inter- 
face. In its crudest form of matching the rotor reactions on a rigid 
support are first determined and then applied to the flexible airframe 
to establish its response. The design goal is to avoid airframe 
resonances to the dynamic rotor inputs. Fig. 2 shotra the results of 
such an airframe response analysis taken from reference 4. The graph 
gives the vertical acceleration amplitude at the pilot station from 
+ 10,000 inch pounds pitching moment input at the rotor hub vs. the 
excitation frequency for various ballast distributions of the AH-56A(AMCS) 
helicopter. While the trend is probably true, the results do not 
include the effects of the ballast changes on the rotor moment input. 

These effects may be substantial. 

A better approximation is obtained when an equivalent hub mass is 
added that is representative of the rotor impedance at the excitation 
frequency. As outlined for example in reference 5, this method of 
establishing the airframe response to rotor dynamic Inputs has its 
limitations. The equivalent masses are different for horizontal and 
vertical motions, and the equivalent mass concept also does not consider 
aerodynamic damping effects. The correct way of dynamically matching 
rotor and airframe is to determine the impedance or mobility matrices 
of both substructures at the interface and to write the compatibility 
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relations for the Interface forces and deflections. Reference 5 
presents a simple example for this method. Since the cruder method of 
exciting the airframe with separately determined rotor Input forces Is 
suggested In textbooks (references 6 and 7) and Is also practiced In 
Industry, It was believed worthwhile to apply the method of Impedance 
matching to a more complex example and to demonstrate the large changes 
In rotor Input to the airframe caused by dynamic rotor/alrframc 
coupling. 

The principle of rotor/alrframe Impedance matching Is seen In 
Fig. 3. A rotor on a rigid support Is schematically shown In Fig. 3a. 
The force amplitude on the support Is 1 q. This symbol should be 
Interpreted as the column of complex valued force and moment ampli- 
tudes that the rotor transmits to the rigid support. Fig. 3b shows a 
rotor supported by an airframe. Here f represents the dynamic force 
column, z Che deflection column Including angular deflections. The 
force f is equal to f^ minus the rotor reaction force due to the 
deflection z at the rotor/airframe Interface. 

f = fg ■ ^ 

where Q„(X) is the rotor Impedance. The same force f acts on the 
R 

airframe or rotor support and produces the same deflection z 

f = 2 (2) 

where Q«(X) is the airframe or support impedance. By inserting Eq. (2) 
Into Eq. (1) one obtains 

f - (I + Q^(X) Qg(X)‘^)'^ fQ (3) 


I Is Che Identity matrix. 
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Rotor and Support Representation by 3 Masses and 2 Springs 

Before deriving the Impedance matrix for an actual rotor, a 
quantitative interpretation of Eq. (3) will be given for the simple 
system defined in Fig. 4. The rotor is represented by two vertically 
moving equal masses m^ Interconnected by a spring with stiffness 
The rotor support is represented by a mass mg and by a spring with 
stiffness Kg. A dynamic force f^ with imaginary frequency X acts on 
the support mass mg if this mass is rigidly held. The impedances are 


now polynomials of X and Eq. (3) can be written 

f/fp = (1 + Q^(X) Qg(X)"^)’^ (4) 

Introducing 

4 “ ^/“r » “s “ ^5) 

it is easy to show that 
Qr(X) = 

Qg(X) = mg(X^ + o)g) (7) 

Inserting Eqs. (6) and (7) into Eq. (4) 

f/fQ * [1 + (mj^/mg)(X^ + u)^/(l + + tOg)]“^ (8) 


Selecting a time scale such that u)_ = 1 we have for f/f « <» the 

K O 

characteristic equation 

X^(l + nij^/mg) + (u)g + 1 + 2 m^/nig) + Wg = 0 (9) 

Since the system Is conservative the roots of Eq. (9) are purely 
Imaginary, X = io). For o)^ = * 0 one has the natural frequency 
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- 0)g/(l + 

For « K^ o 00 the natural frequency Is 

Wjj - (Ug/d + 2 nij^/mg)^^^ (11) 


Eqs. (10) and (11) follow from Eq. (8) for f/^Q “ “• They are also 
directly evident from Fig. 4. In the first case the spring Kg carries 
the mass m^ + m^ and Is uncoupled from the upper mass m^. In the 
second case the spring Kg carries the mass 2 m^ + mg. 


For (0 


R 


1 and 


ni„ 


2.5 the roots of Eq. (9) are only dependent 
on the support uncoupled frequency o)g. They are shown vs. Ug In 
Fig. 5. The higher of the two natural frequencies Is for large oi. 


asymptotlc to the value of Eq. (10) for aj„ * 0. The lower natural 
frequency Is as 3 nnptotlc to the line ^ 1. This asymptotic behavior 

for large cOg or Kg Is also directly evident from Fig. 4. Inserting 
X = iu) into Eq. (8) and introducing ® » “5 “ Wq/w one obtains 

( 12 ) 


f/fg = [1 - - T))/(“s “ 


Vertical Force Excitation 


For a numerical example we select the case of Fig. 5 with ° 1, 

K. 

m /m„ “ 2.5 and stipulate that the forcing function frequency be 
R b 

0 ) = 1.89. Now 0 ) = 1/1.89 = .53 and Eq. (12) becomes 
R. 

f/fg = (1 - 1.53/(mg - 1))"^ (13) 

For oig * 1 when the excitation frequency ui is equal to the uncoupled 
support frequency Ug, the ratio f/fg is zero. Absolute values from 
Eq. (13) are plotted in Fig. 6. It is seen that the input to the 
support mass m„ depends critically on the support stiffness as ex- 
pressed by the ratio = iUg/u. 
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Eq. (13) can also be derived In a different way making use of the 

value a)_„ associated with the coincidence of excitation frequency ui with 
SN 

the coupled system natural frequency At this frequency “ “ 

or f^ = 0 and we have from Eq. (4) 

Qr(u)j,) “ - QgCwjj) (1^) 

Thus £q. (4) becomes 

f/f^ -= (1 - Qg(ujjj) QgCo))”^)"^ (15) 

Inserting the support Impedances from Eq. (7) 

f/fp = (1 - (o)gjj - <4)/(wg - 0)^))'^ (16) 

Stipulating = w and using as before ® 

-2 _2 

f/f^ - (1 - (ujgjj - l)/(o)g - D) ^ (17) 

For the example of Fig. 6 we have “ 3/1-89 ® 1.59 and Eq. (17) 
reduces to Eq. (13). Eqs. (15) and (17) are particularly convenient 
sxace they use only support data and no rotor data, except that the 
support stiffness resulting in a coupled natural frequency of = oj 
muL,L be known. Eq. (17) will be used later in connection with a more 
complete rotor analytical model. 

It is also of interest to deteraine the springforce amplitude f^ 
that the support spring transmits to the base, see Fig. 4. One can 
easily derive the expression 

The absolute \^alues are plotted in Fig. 7 vs. Ug assuming as before 
(i)gj^ = 1.59. The spring force amplitude is smaller than the rotor 
excitation force f when to < 1,1. Fig. 7 also shows the base input 

O o 

4 ^ 
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f /f for f/f =1.0 according to 
SO o 

For low support stiffness this base input is much larger than the 
actual value. Indicating again the error encountered vdien ignoring 
rotor-support coupling. 

For a DAVI type of support (D3mamic Antiresonance Vibration 

Isolator, see for example references 8 and 9) with a large ratio of 

2 

1/a, see Fig. 8, and with an equivalent absorber mass m » m (1/a) 

A 

the force amplitude transferred to the base is 
ihis base force is zero for 

<0g = (m^/(m^ + (21) 

The ratio f /f for m /m^ * 1 and for the same * 1*59 as before is 
s o A S SN 

also shown in Fig, 7. Zero base input occurs at = .7. At this low 
support stiffness the base input without DAVI is already rather low, 
namely f/f - .3. 

From the preceding example it is evident that the excitation of 
the airframe with the rotor force acting on a rigid support can lead 
to large errors in the response. When the excitation frequency coin- 
cides with an uncoupled natural airframe frequency, one would erro- 
neously obtain a large resonance response, while actxially the response 
is modest and no force is transmitted to the airframe. 


Roiling and Pitching Impedance Analysis for a Hingeless Rotor 

The following rotor impedance analysis neglects chordwise blade 
vibrations and is limited to zero advance ratio. The rotor impedance 
is not expected to vary drastically with advance ratio as long as the 
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advance ratio is moderate. Thus the use of rotor Inputs to a rigid 
support at a given advance ratio together with the impedance at zero 
advance ratio should result in a reasonable approximation for the 
vibration prediction. Within the advance ratio range in which a 
constant coefficient multiblade rotor representation is feasible 
(reference 10) it is easy to apply corrections from non-zero advance 
ratio. Only when periodic coefficients are required also in multi- 
blade coordinates will it be necessary to reformulate the rotor imped- 
ance problem* 

The impedance analysis uses blade element transfer matrices 
similar to those described in reference 11. Each blade element has 
constant bending stiffness and is free of distributed loads. The 
stiffness can vary stepwise from element to element. Point loads are 
assumed to act at the element boundaries. Except for the stepwise 
stiffness variation and except for the point loads at the element 
boundaries the analysis? is exact with no further approximations. This 
has the advantage that large beam curvature can be admitted for example 
in the cantilever blade root section without the need of a large 
number of short elements as they are required for the C 81 computer 
program. In reference 11 aerodynamic loads are omitted, so that only 
vacuum modes and natural frequencies are obtained. They are not 
suitable for computing blade or rotor impedances. In reference 11 the 
purely imaginary eigenvalues of the blade characteristic equation are 
found by trial and error, a method not feasible for non-conservative 
systems with complex valued natural frequencies. 

The four steps to extend reference 11 method to our problem are 
the following: 

1. Aerodynamic terms are included in the blade transfer matrix as 
point loads at the element boundaries using quasi steady theory. 

2. Relations between blade root state variables are derived in- 
cluding a dynamic blade pitch term. 

3. Single blade polynomials are transformed into multiblade matrices 
applicable to 3 or more blades per rotor. 

Total system equations are derived by mobility or impedance 
matching at the rotor /airframe interface. 


4 . 
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Blade Elemont Transfer Matrix 


The computation procedure will be outlined here, tfore details 

are given in reference 12. As Indicated in the Nomenclature all 

quantities are non-dlmenslonalized by using the length unit R, the 

time unit 1/n and the mass unit p R, where p is the blade mass per 

o o 2 2 

unit length at the blade root. The force unit then is p.R n , the 
3 2 

moment unit is p^R ^ , etc. A blade element with its inputs at the 
boundaries is shown in Fig. 9. There are 5 state variables; shear 
force S, bending moment M, slope y’» deflection y, and pitch angle 9, 
forming the state vector 


X 



( 22 ) 


The transfer matrix from one element to the next is defined by 


= [E] [F] (23) 

The matrix [F] gives the change in state vector from the right side of 

the mass m to its left side. It expresses the effects of vertical 
^ 2 

inertia force (-^ m^) , of the aerodynamic forces from damping and 

from dynamic blade pitch angle distributed as point forces at the 
element boundaries (-X c^, g^^), and the centrifugal force element from 
the mass m^ (- m. r.). The matrix [F] also expresses the continuity 
relations between the right and left sides of the mass m^^. 

v[ = y;. = y^. ©1 - ©r 


The matrix [F] has the form 
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[F] 


10 0 •!• Xc^) 

0 10 - nij^ r j 0 

0 0 1 0 0 


0 0 0 1 0 

0 f' 0 0 1 


(25) 


The matrix [E] gives the change in state vector from the left side of 
the mass m^ to the right side of the mass m^ 1. It represents the 
exact solution for the bending of a uniform beam under tension. 


1 

0 

0 

0 

0 

^i 

1 

0 

0 

0 

^31 

E32 

^33 

^34 

0 

^41 

^42 

"43 

^44 

0 

0 

0 

0 

0 

1 


where 

E31 * - E^2 “ " (cosh(pap - 1) 

E32 “ “ Pi sinh(p^l^) 

^33 ■ ^44 ■ 

^34 “ ^i si^^^Pi^i) 

E^^ « (1/T^) (1/p^) sinh(p^l^) 

E ^3 = - (1/p^) sinh(p^l^) 
p. - (T./q EI.)^^^ 


(26) 


(27) 


In comparing these expressions with those given in reference 11 one 

must note that X is defined differently leading to the opposite sign 
2 

of X in Eq. (25). Also, the unit of circular frequency in reference 11 
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4 1/2 

is (EIq/PqR ) Here it is the angular rotor speed Q. Introducing 

the transfer matrix 


[El [F] » [<»] 


(28) 


Eq. (23) can be written in the form 


’'iM ■ >♦!! *1 


(29) 


For N point masses one obtains by successive transfer matrix multiplica^ 
tion» beginning at the blade tip, for the blade root quantities the 
following relations 




*11 *12 *13 *14 


"Sl“ 


"*15 



*21 *22 *23 *24 


«1 


*25 


- 

*31 *32 *33 *34 




*35 

^N+1 


*41 *42 *43 *44 


_^1_ 


_^5_ 


(30) 


The blade state variables with subscript 1 refer to the blade tip, 
those with the subscript N+1 refer to the blade root assumed to be at 
the rotor center. The blade pitch angle G is assumed to be constant 
over the blade radius so that the system of Eq. (30) has only 4 equa- 
tions rather than the 5 of Eq. (29). All elements of [<{>] are poly- 
nomials of X. Eq. (30) completes step 1 of the impedance analysis. 

The computations of the aerodynamic coefficients c^ and are straight 
forward and not presented here* For a uniform blade they are propor- 
tional to the Lock number y which in our non-dimensional units is 
given by Y * 3 P a c. 


Blade Root Moment Impedance 


We assume zero vertical motion at the rotor center 


^N+l 


0 


(31) 
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The blade tip bending moment and shearforce is zero 


- 0 


( 32 ) 


The two remaining tip variables, and will be expressed in terms 
of the blade root variables by inserting Eqs. (31) and (32) into Eq. 
(30). 


®N+1 

ss 

•^“13 

'•’14 


y’l 

+ 

*15 

e 

(33) 

*Vl 
— ^ 


•^23 

*24_ 


_yij 


_*25_ 



'vr 

! 

* 

***33 

*34 

1 

1 

y'l 

+ 

*35 

© 

(34) 

0 

i 

'^3 

*44_ 




_*45_ 

1 

1 



Computing 


L^U 


from Eq. (34) and substitution into Eq. (33) one obtains 


after some manipulations for the second of Eq. (33) 
%fl * ^'n+I 


where 


(35) 


P„(X) 



*33 

*34 


*23 

*34 

11 



, P (X) « 




'*43 

*44 

(X 

*43 

*44 


Pg (X) 


23 

-e- 

*43 

'24 

*24 

*44 

25 

*35 

*45 
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I£ the dynamic blade pitch angle 6 is linearly related to the root 
slope, Eq. (35) can be used to determine the blade root bending imped- 
ance ^ pitch-flap coupling expressed by 0 - k 

the root bending impedance is 

^ 06 ) 

Eq. (36) completes step 2 of the impedance analysis. 


Transformation to Multiblade Impedance 

The relations between single blade and multiblade coordinates 
are, when only cyclic terms are retained, see reference 10, 

*k + “ii *k 
\ * ®II \ '*'k 

<“»l>k * \ * “ii »k 

a^, M^, refer respectively to nose down hub tilting angle, nose 
down rotor moment, and nose down cyclic pitch input, 0^^ 

refer to left hub tilting angle, left rotor moment and left cyclic 
input. Eqs. 37 are easily inverted, for example 

b 

- (2/b) 2 (y'N+i'k *k 
k«l 

Eqs. (37) are valid for 3 or more blades, b > 3. 

By differentiating the relation between a single blade variable z 
and its multiblade counterpart z^, z^^ 

z ■ Zj cos t + Zjj sin t (39) 

• .. 2 

with respect to the time t and by replacing z, z by Xz, X z, and z_, 

••2 ^ 

Zj by Az , X z , etc., one can prove a general rule for transforming 
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single blade polynomials into multiblade poljmomlals. According to 
this rule the single blade expression 


P(X) z (40) 

with P(X) a polynomial, is transformed Into the multiblade expression 


■ 

Pj(X)1 Tzj' 

(41) 

- Pj(X) 

Pr(^)J L^ii- 


where Pj^(X) and Pj(X) are 

obtained from 



P(X + i) = Pj^(X) + 1 Pj(X) (42) 


by separating real from Imaginary terms. To apply this rule we write 
Eq. (35) as 

Vl “ ■ y'N+1 ^ + 1) W3) 

Using the definitions of Eq. (37) the multiblade equations are 








\ ' 



-^II- 




-“ll- 


(44) 




1 

o 

M 

1 



L^xiJ 


In the absence of cyclic pitch inputs from the control system we have 
according to Fig. 10 


01 = 


■'ll 


= X 


II 


(45) 
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Insertlng Eq. (45) into Eq. (44) one obtains the rotor Impedance 
matrix in roll and pitch by premultiplying Eq. (44) with 




( 46 ) 


This completes step 3 of the impedance analysis. 


Coupled Rotor/Support Analysis 

A coupled rotor/support analysis is performed for the system shown 
in Fig. 11. The rotorshaft is assumed to be rigid and connected to a 
rigid housing that is supported by an up-focusing mount with focus on 
the rotor center. Thus the hub is rigidly restrained against hori- 
zontal and vertical motions. The elastic restraints in pitch and roll 
are indicated by the horizontal springs attached to the shaft housing. 
There also is a gravitational restraint (pendulum effect) that is to 
be included in the support stiffness and in the support natural fre- 
quency. 

3 2 

The unit for the moments M_, is p R . With the time unit 

1 II o 

1/fi one then obtains the rotor pitching and rolling moments 

(b/2) Pq Mj , (b/2) p^ R^ Mjj (47) 

For uniform blades 


p R^ = 3 I, (48) 

O D 

The rotor support d3mamic equations without support damping are 


2 2^ 


— > — 


P — T 

(u - X 0 




“l 

I 





^ 2 ,2 



= 3(b/2) (Ig/I) 


0 0)-- — X 





II 


II 
L - 


II 


( 49 ) 
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6y Inserting 


ftj 


from Eq. (49) into Eq, (44) one obtains with Eq. (45) 


a set of homogeneous equations for a 


that are of the form 



(50) 


The coupled system characteristic equation is then given by 


= 0 (51) 

The mode shapes are computed by inserting intc Eq. (50) one of the 
eigenvalues obtained from Eq. (51). One then finds that all modes are 
either regressing or progressing in the rotating reference system. 

This completes step 4 of the analysis as far as coupled system natural 
modes and frequencies are concerned. The ratio of the rotor moment on 
the flexible support over that on a rigid support will be determined 
from Eq. (17) as a function of the support uncoupled natural frequency. 
\11 that is required for this purpose is the knowledge of the uncoupled 
support frequency for which the excitation frequency coincides 
with the coupled system natural frequency. 



Computational Limitations 


The computations were performed on the IBM-360/65 computer using 
double precision (16 digits). Single blade computations were made for 
5, 8, 10 15 and 20 elements per blade. The number of 20 elements was 

found to be too high for the 16 digit precision used. With 10 blade 
elements no computational difficulties were encountered, provided that 
the evaluation of the 3 by 3 determinant in Eq. (35) was numerically 
optimized by writing 
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I...1 “ '^33^^24^5 “ ^4'<’25^ *34^‘*’23^5 ‘ '**43'^25^ " ‘<’35 ^'**23^4 " ‘<’43^24^ 

(52) 

whereby the set <l> 32 > <t> 3 ^ represents the column with the highest 

values of the elements. A convenient check for adequate computer 
precision consists of looking at the coefficients of 
etc., in P^, whereby nmax is the highest power that should occur 
theoretically in this polynomial. With N point masses nmax « 2N-2 for 
Pg. If the computer precision is adequate, the coefficients of 

^nmax+1^ ^nmax+2^ etc., are several orders smaller than those of 
nni£ix 

X • The accuracy of the computation depends on how the determinant 
is evaluated, the method of Eq. (52) giving the best results* 

Once the single blade problem could be solved without difficulties, 
the coupled system solutions posed no further obstacles. The computer 
accuracy depends on the highest eigenvalue considered and does not 
suffer when the number of eigenvalues is approximately doubled as for 
the coupled system, as long as the highest eigenvalue remains approxi- 
mately the same. The results presented here were obtained with 8 
masses per blade. As compared to an analysis with 10 elements per 
blade the first 3 eigenvalues shown here were found to have less than 
1% error. The equivalent CPU time to obtain a complete set of eigen- 
values for the coupled system was 4 seconds. For 55 cases computed in 
one run the CPU time was 34 seconds. 

The introduction of more blade degrees of freedom - in-plane 
bending and torsion - and of a much higher order of the airframe 
impedance polynomials may not lead to computational difficulties, if 
the highest eigenvalue of the system is not substantially increased. 

If numerical difficulties do occur, the method of reference 13 would 
be helpful. Accord.lng to this method the order of the characteristic 
equation is reduced by modalizing and truncating the component mobility 
matrices and by applying a correction term to the truncated expressions. 
The method has been successfully applied to the problem of Eq. (51). 

The results will be reported elsewhere. 
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Coupled Rotor and Support Modes 

The two parameters that determine the dynamics of the uniform 

blade are the non-dlmenslonal blade bending stiffness q and the Lock 

1/2 

number y. The centrifugal tension force occurring In p = (T^^/q) 
of Eq. (27) can be easily computed for each blade station 1. The 

blade Lock number occurs as factor In the aerodynamic terms c and g 

- 1/2 ^ ^ 
of Eq. (25). Since q Is proportional to the rotor angular speed 

It can be used as non-dlmenslonal speed parameter. We select for the 

- 1/2 

numerical example q =18 and y = 5. This results In a blade 
cantilever first natural frequency In vacuum of 1.06, which Is a 
realistic value for a hingeless rotor. For the first 3 cantilever 
blade bending modes one then obtains the following complex eigenvalues, 
whereby R and P denote regressing and progressing modes In I'he rotating 
system respectively. 

- 1/2 

Rotor Alone, q = 18, y « 5 
Natural Frequency 5 ± lui 


Flap-Bending 

Hingeless Rotor 

Hinged Rotor 

Mode 

5 1 

0 ) 

C 

0 ) 

IR 

1 

0 1 

0 


.05 




-.311 


IP 

-.616 

1.96 


1.95 

2R 

-.327 

1.57 


1.57 




-.257 


2P 

-.193 

3.62 


3.57 

3R 

-.240 

3.87 


3.80 




-.229 


3P 

-.214 

5.87 


5.80 


For comparlslon purpose the eigenvalues for the hinged rotor have been 
added, assuming hinges at the rotor center. Since there is no coupling 
between blades it is sufficient to perform a single blade analysis. 
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whereby from Eq. (36) (X) * 0. This gives the imaginary parts in 

the rotating system. Those in the non-rotating system are obtained by 
0 ) “ 01 ^ ± 1. The imaginary parts of the eigenvalues are almost the 
same for hinged and hingeless rotor. The interblade coupling for the 
hingeless rotor affects only the real parts, that is the damping of 
the modes. While this damping is the same for progressing and regressing 
modes of the hinged rotor, it differs for the hingeless rotor. The 
largest difference is for the first progressing mode where the damping 
of the hingeless rotor is twice that of the hinged rotor. 

The hinged rotor does not couple with the support shown in Fig. 11 . 
The coupling of the hingeless rotor is determined by the support 
natural frequencies and by the parameter 3(b/2) I^/l of Eq. (49). 

We select the blade number b » 4 and the ratio of blade moment of 
inertia over support moment of inertia of Ijj/I * 5. We further stipulate 
(D^ = = o)g which means equal support stiffness in pitch and in 

roll. When performing the coupled system analysis it turns out that 
only the two hingeless rotor modes IP and 2R substantially couple with 
the support modes, while the higher rotor modes 2P, 3R, 3P are not 
much affected by the coupling with the support. The coupling of the 
IR mode with the support leaves the imaginary part almost at zero but 
introduces substantial damping, resulting In a almost aperiodic mode. 

When varying the support stiffness expressed in uig and keeping 
all other parameters (q, y* b, Ij^/I) constant, one obtains for the 
coupled IP modes the 5 and u) values shown in Fig. 12. The uncoupled 
values are indicated in dash lines. There is a region of instability 
indicated by a positive real part 5 between an uncoupled support 
frequency of 1.5 and 3.0. Soft mounting a hingeless rotor can thus be 
dangerous. Note that we are not considering so called air resonance 
where the In-plane regressing mode becomes unstable but that we have 
here a hingeless rotor progressing llap mode instability on its soft 
support, that is usually called whirl flutter. Similar results would 
be obtained for a hinged rotor with sufficiently large off-set of the 
hinges from the rotor center. 

The t, and values of the coupled 2R mode are shown in Fig. 13. 

This graph is quite similar to Fig. 12 for the IP mode, except that 
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there Is no unstable region of support stiffness. The damping ratio 
of the coupled support modes, whether positive or negative, approxi- 
mately given by c/u. Is not large. This justifies the neglect of 
damping In the next section. 


Hingeless Rotor Moment Excitation 

In order to find out the effect of rotor and support coupling on 
the moment transferred from the rotor to the support, we apply Eq. C17) 
to the pitching and rolling moment. 

M/M^ - [1 - - l)/((aig/(u)^ - 1)]“^ (53) 

For a A bladed rotor the excitation frequency is oo ® 4. For both the 
IP mode (Fig. 12) and the 2R mode (Fig. 13) the uncoupled support 
natural frequency u) for which the excitation is in resonance with 

wN 

the coupled rotor-support natural frequency is 





“SN ■ 

(54) 

Thus from Eq. 

(53) 





M/M “ 
o 

(1 + 

• lO/CCujA)^ - 1))"^ 

(55) 

with 

M/M = 
o 

0 

for ojg = 4.0 

(56) 


M/M - 
o 

00 

for aig = 3.8 

(57) 

The relation 

of 1 M/M^ I vs . 

■ “s 

is shown in Fig. 14. Note 

the very 


steep descent of the curve between * 3.8 and * 4.0. Same as for 
the example of vertical vibrations shown in Fig. 6, the value = 1 

is reached asymptotically for large values of (Og representing high 
support stiffness. Zero moment transfer from the rotoi to the support 
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is obtained for Wg ® 4*0 when the uncoupled support natural frequency 
is in resonance with the excitation frequency. If damping had been 
considered, would not have reached a value of zero at ojg = 4 

but rather a minimum value different from zero. 

The spring moment transferred to the base is according to Eq. (18) 

M /M * (58) 

so S S SN 

This ratio is plotted in Fig. 15 vs. Wg and shows that it is less than 
one for a>g < 2.7 where the system is unstable (Fig. 12). Using DAVI's 
the base moment could be made almost zero for an uncoupled support 
frequency of < 3.8. This condition could be unstable. 


Conclusion 


1. The method of rotor induced vibration prediction by applying the 
rotor d 3 mamic forces and moments acting on a rigid support to the 
flexible airframe can lead to large errors of either over or 
under prediction of vibrations. This method also cannot account 
for possible dynamic instabilities from rotor-support coupling. 

2. The rotor dynamic forces and moments acting on a rigid support 
must be corrected by a term that has as a factor the rotor imped- 
ance taken at the rotor/airframe interface. 

3. A practical way of determining the hingeless rotor impedance 
matrix for the pitching and rolling moments has been developed by 
first performing a finite element single blade analysis including 
aerod 3 mamic terms and by then applying a simple transformation to 
multiblade impedances. 

4. Extensions of the rotor impedance analysis are desirable with 
respect to the following items; 

4.1 From zero to moderate advance ratio by using constant 
coefficient multiblade equations. 

4.2 From moderate to large advance ratio where periodic coef- 
ficients in multiblade equations are necessary. 



24- 


4.3 From a 2 by 2 impedance matrix for roll and pitch to a 
higher order impedance matrix to include linear hub forces 
and deflections. 

4.4 From rigid to flexible blades in chordwise bending and 
torsion » using the transfer matrices developed in reference 
14 and extending them to include aerodynamics. 

4.5 From a simple rotor support to a complete airframe. 
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